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Abstract 



We present a solution to the cosmological problem encountered in (su- 
persymmetric) grand unified theories due to copious monopole production 
at the end of hybrid inflation. By employing thermal inflation "driven" 
by the U(l) axion symmetry, the superheavy monopole flux can be nat- 
urally suppressed to values that should be accessible to dedicated large 
scale experiments. The U(l) axion symmetry also helps generate the right 
magnitude for the /i term of the minimal supersymmetric standard model. 
An important by-product is the predicted existence of stable or very long- 
living fermions possessing intermediate scale masses of order 10 12 GeV. 
Their presence is required for implementing thermal inflation, and their 
stability is due to a Z2 symmetry. They may constitute a sizable frac- 
tion of cold dark matter, and possibly help explain the ultra-high energy 
cosmic ray events. The rest of cold dark matter may consist of axions. 
Although our discussion is carried out within the framework of supersym- 
metric SU(4) C x SU(2)l x SU(2)r , it can be extended to other grand 
unified gauge groups such as SU(3) C x SU(3) L x SU(3) R or 50(10). 
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The great advantage of hybrid inflation [|1] is that, in contrast to previous inflationary 
schemes, it can reproduce the observed temperature fluctuations of the cosmic microwave 
background radiation with natural values of the relevant coupling constant. Moreover, 
this inflationary scenario is almost automatically realized [0,0] in supersymmetric (SUSY) 
grand unified theories (GUTs). However, in trying to apply it to GUTs which predict the 
existence of magnetic monopoles, a cosmological disaster is encountered. Hybrid inflation 
is terminated abruptly when the system reaches an instability point on the inflationary 
trajectory and is followed by a 'waterfall' regime during which the spontaneous breaking 
of the GUT gauge symmetry takes place. The appropriate Higgs fields develop their 
vacuum expectation values (vevs) starting from zero and they can end up at any point of 
the vacuum manifold with equal probability. As a consequence, monopoles are copiously 
produced |4| by the Kibble mechanism leading to a cosmological catastrophe. 

Possible solutions to this monopole problem have been proposed [|,|J. They rely on 
introducing the leading non-renormalizable term in the standard superpotential |2|] for 
hybrid inflation. (For different resolutions of the problem see Ref. |7|].) In Ref. Q, the tri- 
linear coupling of this standard superpotential was eliminated by a discrete symmetry and 
was replaced by the leading non-renormalizable term. The system, from the beginning 
of inflation, follows a particular valley and ends up at a particular point of the vacuum 
manifold. Thus, no monopoles can be produced. The inflationary trajectory possesses 
a classical inclination driving the inflaton towards the SUSY vacua and the termination 
of inflation is smooth. In Ref. ||, both the trilinear and the leading non-renormalizable 
couplings were kept revealing a quite different picture. The inflationary trajectory is 
classically flat and, thus, radiative corrections are needed for driving the inflaton. 
The termination of inflation is abrupt. Nevertheless, there is no monopole production 
since the GUT gauge symmetry is already broken during inflation. Both models predict 
complete absence of monopoles which may be disappointing for the experimenters. 

In this letter, we propose an alternative solution to the monopole problem of hybrid 
inflation which may yield a measurable monopole flux in our galaxy. The idea is to keep 
the original SUSY hybrid inflationary scenario unaltered and try to dilute the monopoles 
by invoking a subsequent thermal inflation associated with an intermediate mass 
scale. (The main ideas underlying thermal inflation have been presented in Ref. H. The 
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term thermal inflation was coined in Ref. ||, which further elaborated on the scheme.) 
It is then natural to identify this scale with the one at which the Peccei-Quinn (PQ) 
symmetry [ID is broken and also use [[TTJ it for generating the fi term of the minimal 
supersymmetric standard model (MSSM). Although of much wider applicability, our 



mechanism is presented in the context of the SUSY Pati-Salam (PS) model |[L2| , which 
is one the simplest unified schemes possessing Jl3| monopoles. 

Our mechanism has an interesting by-product. In trying to make thermal inflation 
possible, we are led to the introduction of a number of superfields coupled to the field 
which breaks spontaneously the PQ symmetry. These fields possess intermediate scale 
masses, with the lightest ones being either stable or very long-living as a consequence of 
a Z2 symmetry. Their fermionic components may constitute a sizable fraction of the cold 



dark matter in the universe, and possibly help explain [14] the ultra-high energy cosmic 



rays |15[ . Axions, of course, also may contribute to the cold dark matter fraction. 

We consider the SUSY PS model with gauge group G PS = SU(A) C x SU(2) L x 
SU(2)r . The left-handed quark and lepton superfields are accommodated in the repre- 
sentations Fi = (4, 2, 1), F? = (4, 1, 2), where i = 1,2,3 is the generation index. The two 
electroweak Higgs superfields belong to the representation h = (1,2,2). The PS gauge 
group can be spontaneously broken to the standard model gauge group by a conjugate 
pair of Higgs superfields H c = (4, 1, 2), H c = (4, 1, 2) acquiring non- vanishing vevs along 
their right-handed neutrino directions. This can be achieved by introducing a gauge sin- 
glet superfield 5" with two (renormalizable) superpotential terms: a term linear in S and 
a trilinear coupling of S to H c , H c . The resulting scalar potential automatically pos- 
sesses an in-built (classically) flat direction along which hybrid inflation can take place 
with the system driven by an inclination from one- loop radiative corrections ||. Gps 
is restored along the inflationary trajectory and breaks spontaneously only at the end of 
inflation when the system falls towards the SUSY vacua. This transition leads H to a 



cosmologically catastrophic copious production of doubly charged monopoles ||13|| . The 
monopoles could be diluted to an acceptable level if the primordial hybrid inflation is 
followed by thermal inflation [|8],|9[]. This inflation, associated with an intermediate mass 
scale, is terminated at cosmic temperatures of the order of the electroweak scale and 
generates only a moderate number of e-foldings. 
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Thermal inflation could be "driven" by the PQ symmetry which is spontaneously 
broken at an intermediate scale. Moreover, in Ref. we showed that the PQ sym- 
metry, which solves the strong CP problem, can also be used to generate the /i term 
of MSSM with the desired magnitude. More specifically, we introduced a pair of gauge 
singlet superfields N, N with non-zero PQ charges and a (non-renormalizable) super- 
potential coupling N 2 N 2 . The resulting scalar potential of these fields, after soft SUSY 
breaking, possesses a non-trivial minimum which, under certain circumstances, is the 
absolute minimum, with N, N acquiring intermediate scale vevs, thereby breaking the 
PQ symmetry. The fi term is then generated via the superpotential coupling N 2 h 2 . The 
trivial extremum (at iV = N = 0) turns out to be a local minimum separated from the 
PQ minimum by a sizable potential barrier. This situation persists at all cosmic temper- 
atures after the primordial reheating which follows hybrid inflation, as has been shown 
by including the one-loop temperature corrections to the potential. Thus, a successful 
transition from the trivial to the PQ minimum cannot be realized. We had to assume 
that the system, after the primordial reheating, already emerges in the PQ vacuum. In 
other words, there is neither PQ transition nor thermal inflation in this case. 

In order to make thermal inflation possible, we must turn the trivial local minimum 
of the zero-temperature potential into a saddle point. More specifically, the positive soft 
mass 2 term of N should become negative. This can be achieved radiatively and requires 
strong couplings of iV to a number of superfields. To this end, we introduce n superfields 
D a (a = 1, 2,...,n) belonging to the representation (6,1,1) of Gps with superpotential 
couplings NDD. However, these color (anti)triplets acquire intermediate scale masses 
after the PQ breaking, which could prevent the unification of the MSSM gauge coupling 
constants. To restore gauge unification, we include an equal number of superfields H a 
belonging to the (1, 2, 2) representation with superpotential couplings NHH . Note that 
the negative mass 2 term of N, which is successfully generated by invoking these extra 
superfields, gives rise to an additional problem. All the higher order terms in the scalar 
potential, which are derived from N 2 N 2 after soft SUSY breaking, involve both N, N. 
Thus, these terms vanish along the N direction, and the potential becomes unbounded 
below due to the negative mass 2 term of N. Fortunately, there is a simple way out from 
this "runaway" problem. Replacing N 2 N 2 by N 3 N, we generate a |iV| 6 term in the 
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potential which prevents the runaway behavior along the N direction. 

We still need to include some extra couplings and superfields to obtain a phenomeno- 
logically viable scheme. In particular, we must introduce quartic superpotential couplings 
of H° to Ff. These couplings generate intermediate scale masses for the right-handed 
neutrinos and, thus, seesaw masses for the light neutrinos. The inflaton then decays 
into right-handed neutrinos via the same couplings. Finally, in order to give superheavy 



masses to the down quark type components of H c , H c , we include [0] an SU(A) C 6-plet 
superfield G = (6, 1, 1) with superpotential couplings GH C H C , GH C H C . 

The superpotential of the model, which incorporates all the above couplings, is 

N 2 h 2 N 3 N 

W = nS(H c H c - M 2 ) + Ax + A 2 + a a ND a D a + f3 a NH a H a 

mp mp 

+\ lJ F l c F j h + li ^^F c i F<: + aGH c H c + bGH c H c , (1) 

where M is a superheavy scale and mp = Mp/^/8n « 2.44 x 10 18 GeV is the reduced 
Planck mass. Here, we chose a basis in the D a , H a space where the coupling constant 
matrices a and f3 are diagonal. Assuming that, at a more fundamental level, the D's and 
if's originate from SO (10) 10-plets, we can obtain the restriction a a — (3 a (a — 1, 2, n). 
Note that M, k, Ai,2, a a , Pa, a and b can be made positive by field redefinitions. 

In addition to Gp$, the superpotential in Eq.(|l|) possesses two continuous global 
(anomalous) symmetries, namely a R symmetry U(l)p and a PQ symmetry U(1)pq . 
The R and PQ charges of the superfields are assigned as follows: 

R : H c (0), H c (0), S(A), G(4), D(l), H(l), F(2), F c (2), N(2), N(-2), h(0); 
PQ : H c (0), H c (0), S(0), G(0), D(l), H(l), F(-2), F c (0), N(-2), N(Q), h(2). (2) 

Note that the R charge of W is 4. Although it is not necessary, we also impose, for 
simplicity, a discrete Z™ p symmetry ("matter parity"), under which F, F c change sign. 
Additional superpotential terms allowed by the symmetries of the model are 

FFH C H C NN, FFH c H c h 2 , FFH C H C NN, FFH c H c h 2 1 F C F C H C H C , (3) 

modulo arbitrary multiplications by non-negative powers of the combinations H C H C , 
(H c ) 4 , (H C ) A (this applies to the terms in Eq.(|I]) too). Note that, without the Z™ p 
symmetry, the coupling DHFH C would also be present in Eq.(H). 



Instanton and soft SUSY breaking effects explicitly break U(1)r x U (1) pq to a discrete 
subgroup. It is then important to ensure that this subgroup is not spontaneously broken 
by the vevs of N, N since otherwise cosmologically disastrous domain walls will be 
produced in the PQ transition. This requirement implies that the number n of D's and 
if's must be 5 or 7. Moreover, in both these cases, the subgroup of U(1)r x U(1)pq 
left unbroken by instantons and SUSY breaking coincides with the one left unbroken 
by (N), (N), and is a Z 2 x Z 2 generated by (e in/2 ,e in/2 ) and (l,e i7r ). (Note that the 
element (e l7T , e l7T ) is equivalent to the identity element since it leaves unaltered all the 
superfields of the theory.) Combining appropriately these Z 2 s with Z™ v and the Z 2 
center of SU(2)l, we obtain two equivalent Z 2 s under which the D J s or the i?'s change 
sign. It is interesting to note that, even if only N develops a non-zero vev (see below), 
absence of domain walls still implies n = 5 or 7. The soft SUSY breaking terms respect 
the symmetry Z 2 x U(1)pq, where the non-anomalous Z 2 is generated by (e i?r//2 , e" 1 "/ 2 ). 
It is then obvious that further breaking of U(1)pq to Z 2 by the vev of N can solve the 
strong CP problem. Thus, we have the option to keep (N) = 0. The symmetries which 
survive after instanton effects, in this case, are the same as in the (N) ^ case. 

Let us note that baryon and lepton number violations arise from the last three super- 
potential terms in Eq.(|3|), the last two couplings in Eq.(|IJ) and the combinations (i? c ) 4 , 
(H ) 4 , in complete analogy with Ref. ||. The proton is practically stable. 

The part of the tree-level scalar potential which is relevant for the PQ (and R sym- 
metry) breaking can be derived from the superpotential term N 3 N and, after soft SUSY 
breaking mediated by minimal supergravity, is given by (compare with Ref. |TT|| ) 

V PQ = ml /2 (\Nf + |Af + Xj W (|Af + 9|iVT) - \A\X 2 ™^|Af ) , (4) 
1 \ {m 3 / 2 mp) z m 3/2 m P J 

where m 3 / 2 is the gravitino mass and A is the dimensionless coefficient of the soft SUSY 
breaking term which corresponds to the superpotential coupling N 3 N. Here, the phases 
(p, 9 and 9 of A, N and N are taken to satisfy the relation tp+39+9 = n which minimizes 
the potential for given values of \N\, \N\. 

As a consequence of the couplings of iV to the -D's and if's, its mass 2 in Eq.(fD can 
easily turn radiatively to negative values at lower energy scales (see also Ref. |17| ). To 
see this, we consider the one-loop renormalization group equations for the soft masses 
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m N, m D and mil of the scalar fields N, D and H (see Ref. |T 

yL ~d K = 3nYD ( 2m D + + 2nY H (2m 2 H + m 2 N ) , 
^drnl = y D ( 2m 2 + m 2^ = y H ( 2m ^ + (5) 

where /1 is the running energy scale, Yd = a 2 /2n 2 , Yh = /3 2 /2ti 2 , with a = a a , (3 — (3 a 
(a = 1,2, ...,n). Here we take universal "asymptotic" soft scalar masses (m^ = mo = 
mu = ^3/2 at the GUT scale) and assume, for simplicity, that all the a's (/3's) are 
equal and, thus, all the .D's (-ff's) have the same soft mass at all scales. This system of 



equations possesses [18| a non-trivial fixed point given by 



2 2 2 2(57*1 9 / \ 

m N = -2m D = -2m H = — — — m 3/2 , (6) 



5n + 2 



and admits the solution 



m 2 N 2(5n-l) 15n ( ft \( 5n + 2 ) Y 



m 2 /2 5n + 2 + 5n + 2\M G J 

m 2 D m 2 H _ 5n - 1 3 / fx \ ( 5n + 2 ) Y 
m l/2 ~ m l/2 ~ 5ra + 2 5n + 2\M G J 

For simplicity, we further assumed that a = (3 (and, thus, Y D = Y H = Y), and we ignored 
the running of these coupling constants. Taking strong Yukawa couplings a = (3 = 1, 
we see that, already at fi/Ma ~ 10~ 2 , the second terms in the right hand sides of 
the equalities in Eq.(^) are much smaller than 1% of the first terms. Thus, after the 
primordial reheating, the soft masses can be taken equal to their fixed point values. 

The radiatively improved zero-temperature scalar potential is given by Eq.(^j) with 
the |iV| 2 term acquiring an extra negative factor — 2(5n— l)/(5n + 2). The trivial (local) 
minimum of this potential (at N = N = 0) then becomes a saddle point and the absolute 
minimum necessarily lies at a non-zero value of N. It should be noted, however, that N 
is also non-zero at the absolute minimum only if A 7^ 0. To simplify the presentation, 
we choose A = 0. The absolute minimum of the potential then lies at 

., .. f a /10n -2\ 1/4 (m 3/2 m P ) 1/2 
The constant energy density carried by the trivial vacuum is given by 
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Vo . 2 (^lf-'^fL, (9) 

V15n + 6/ A 2 
and is responsible for driving thermal inflation. 

The one-loop temperature corrections to the PQ potential can be calculated by em- 
ploying the formalism of Ref . |19[ . The absolute minimum of the resulting temperature- 



dependent effective potential at high cosmic temperatures T lies at N = N = 0. So, 
after the primordial reheating following hybrid inflation, the system emerges in the triv- 
ial vacuum. The temperature correction to the mass 2 term of the field N is 

5V T « nT 2 (3a 2 + 2[3 2 )\N\ 2 } (10) 

for \N\ <C T. Here, we considered only the main contributions which originate from loops 
with the superfields D or H circulating. As T approaches the critical temperature 



/ 10n-2 

Tc "V(5n 2 + 2n)(3^ + 2^) ™ 3/2 ' (U) 
the potential barrier separating the trivial and PQ vacua becomes vanishingly small and 
the PQ transition takes place close to T c M. For a period preceding this transition, the 
energy density of the trivial vacuum dominates over the radiation energy density and 
the universe undergoes a mild inflationary phase known as thermal inflation ||[|]. This 
is terminated at the PQ transition where the system enters into an oscillatory phase 
about the PQ minimum. The coherently oscillating field N (thermal inflaton) with mass 
m in fi = 2[(5n — l)/(hn + 2)] 1 / 2 m 3 / 2 eventually decays, via the superpotential coupling 
N 2 h 2 , to a pair of Higgsinos thereby reheating the universe. Of course, this decay is 
possible only if the Higgsino mass fi = \\j 2 jmp is smaller than m^/i/2. The decay 
width is T m (2\ff 2 /-nm 2 p ){l — ^) x l 2 m in ^ where e = 2fi/m in fi (0 < e < 1). The maximal 
reheat temperature is achieved at e 2 = 2/3, which maximizes T, and is given by 



(12) 



7rg* {l r )mp 

where g*(T r ) « 89.75 is the effective number of degrees of freedom at the reheat temper- 
ature T r ~ GeV (see below). 

The relative monopole number density um/s (um is the number density of monopoles 
and s the entropy density) remains |20| essentially constant for temperatures between the 
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primordial reheat temperature Tr < 10 12 GeV (see below) and the critical temperature 
T c , where the vacuum energy density Vq is transferred to the oscillating inflaton field. 
The initial number density of thermal inflatons is n in fi ps Vo/m in fi . So, at T c , 



n M n M 2% 2 g if (T c )T^m inf i 



(13) 



n in fi s 45Vo 

where g*(T c ) ps 105.75 + 35n/2, for T c ~ 70 GeV (see below), is the effective number of 
degrees of freedom just before the PQ transition. Note that, for temperatures between 
Tr and T c , the fermionic components of the .D's and H's are massless. The ratio in 
Eq.(|T3"D remains J2T] practically unaltered until T r , where um/s = (nM/ninfi){ninfi/s) 
with rii n fi/ s ~ 3T r /4rrii n fi in the instantaneous inflaton decay approximation. Combining 
this with Eq.(^), one obtains the dilution factor of the relative monopole number density 
during the PQ transition and subsequent reheating: 

f^L\ (T \ ~ (^L\ (T \ * 2 9*{T c )T*T r 



V s J V s J 30l/ 

At the primordial reheating, we have um/s ~ {tim/^infl) (3Tr/ Ami n f l) ; where 
m iNFL = V2kM (uinfl) is O the inflaton mass (number density) for hybrid inflation. 



The ratio um /uj^fl remains |ZT| essentially unaltered after the production of monopoles 
at the end of hybrid inflation until Tr . The initial number density of magnetic monopoles 
can be estimated by the Kibble mechanism || which gives tim ~ (3p/4:7r)m% NFL , 
where p ~ 1/10 is a geometric factor. This, together with uinfl ~ ^7 nfl/ m iN f l , 
where Vinfl — k 2 M a is || the constant energy density driving hybrid inflation, yields 
n M /n INFL ps 3pm 4 INFL /4:7iV INFL and, thus, 

(^pa.t^i. (15) 

The fraction of the inflationary energy density which goes into magnetic monopoles right 
after the end of hybrid inflation is ~ 3pm^ NFL mM /Vinfl — SpuniM/ \f2nM <C 1, for all 
relevant k's. Here tum ~ ^trM/go is the monopole mass, with go being the GUT gauge 
coupling constant. The final um/s can then be found from Eqs. flLl]) and fllBD . 

The primordial inflaton decays into a pair of right-handed neutrino superfields v c 
with mass M u c = rjmiNFL/2 (0 < r/ < 1) via the superpotential coupling H C H C F C F C . 
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T « ~ 1 ;= i /4 ,~ ttttt; » ( 16 ) 



The decay width is « (1/8tt) (Af I/ c/M) 2 (l - r? 2 ) 1/2 mi nfl- Note that M^c also originates 
from the coupling H C H C F C F C and, thus, cannot be bigger than about 2M 2 /mp. The 
primordial reheat temperature is given by 

3 l/2 5 l/4 K l/2 m V2 Mye(1 _ ^ 2)1/4 

2^2vr^ 1/4 (T R )M 1 /2 

where g*(Tp) ~ 234.25 + 75n/2 since, at Tr , the bosonic D's and fTs are still massless. 

For each k, the value of M can be found by simultaneously solving Eqs.(2) and (5) 
of Ref. [^] with their right hand sides divided by an extra factor of 2 due to the fact 
that H c (H c ) contains four SU(2)r doublets. We will take, for defmiteness, k = 4x 10~ 3 
which gives M ~ 9.57 x 10 15 GeV. With these values, hybrid inflation ends only when 
the inflaton field S is infinitesimally close to the instability point (at \S\ = M) of the 
inflationary trajectory, as one can easily deduce from the slow roll conditions. Moreover, 
our present horizon scale crosses outside the inflationary horizon at \S\ ~ 2.63M. From 
Eq . (|i~6D , we find that Tr can take all the values up to about 10 12 GeV. The maximal Tr 
is obtained at rf = 2/3 which maximizes the decay width of the primordial inflaton. The 
corresponding value of M u c turns out to be much smaller than 2M 2 /mp. It is important 
to note that the stringent gravitino constraint [24] on Tr (< 10 9 GeV) is alleviated here 
since the primordial gravitinos suffer considerable dilution during thermal inflation. 

We will now estimate the present D and H particle abundance. Due to the two 
unbroken ^'s under which the -D's and if's change sign independently, these particles 
can only annihilate in pairs. The annihilation processes remain in thermal equilibrium 
at all temperatures between Tr and T c . Moreover, at temperatures higher than T c , 
where the PQ (thermal) transition takes place, the fermionic D's and if 's are massless. 
On the contrary, the scalar components of these fields acquire soft SUSY breaking and 
temperature-dependent masses. Consequently, for m.3/2 ^> T c (see below), the number 
density of these bosons at temperatures close to T c is suppressed, and we thus ignore 
them. The number density of the fermionic D's (H's), for temperatures just above T c , 
is n D (T c ) w 21nC(3)T c 3 /2vr 2 (n H (T c ) w 14nC(3)T c 3 /27r 2 ), where C(3) « 1.2021. 

After the completion of thermal transition, the D (H) superfields acquire inter- 
mediate scale masses \/2af a (\/2/3/ a ). Their total relative contribution to the energy 
density of the universe, immediately following the transition, is given by Qdh{T c ) = 
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Pdh{T c )/ pinfi{T c ) ~ 7n(3a + 2/?)C(3)/ a T 6 3 /v / 27r 2 V < 1, where p infl is the energy density 
in thermal inflatons. VLdh remains essentially constant until T r since pair annihilation of 
the fermionic -D's and if 's is frozen already at T c where these fermions acquire their inter- 
mediate scale masses. Between T r and the equidensity temperature T eq « 2.5 x 10~ 9 GeV 
(for present value of the Hubble constant H rs 65 km/s Mpc), where matter and ra- 
diation have equal energy densities, Qdh is enhanced by a factor T r /T eq and remains 
practically unaltered thereafter. So the present abundance of these particles is 

7n(3a + 2/?)C(3) / a T 3 T r 
V2V ^ ^ ' (17) 

This can easily be of order unity so that the D and H fermions with intermediate scale 
masses constitute a sizable fraction of the cold dark matter in the present universe. 

We are now ready to proceed to a numerical example. We choose the gravitino mass 
m 3/2 — 300 GeV and the number of D's and if's n = 7. The thermal inflaton mass is 
then rriinfi rs 575 GeV and the Higgsino mass \x ~ 235 GeV (recall that we take e 2 = 2/3). 
The critical temperature for the PQ transition can be evaluated from Eq. (|TT| ) and turns 
out to be T c « 69 GeV for coupling constants a = (3 = 1. We further take A2 = 2 x 10~ 3 . 
From Eqs.(P) and (§), we then obtain the axion decay constant f a ~ 7.57 x 10 11 GeV and 
the vacuum energy density which drives thermal inflation Vq ~ 3.16 x 10 28 GeV 4 . The 
parameter Ai is evaluated from the Higgsino mass and comes out to be ~ 10~ 3 . Eq. (|I~2]) 
then yields T r m 2.81 GeV for the reheat temperature after thermal inflation. 

The present relative monopole number density is estimated from Egs-flU]) and flT5]). 
It turns out to be tim/s ~ 4.6 x 10~ 41 Tr for the chosen value of k (= 4 x 10 -3 ). This 
implies that Trs of order 10 10 GeV, which can be naturally obtained from Eq.(|l6|) by 
appropriately adjusting M u c (or 77), lead to um/s ~ 10 -30 . This corresponds to the 
well-known Parker bound for the monopole flux in our galaxy derived from galactic 
magnetic field considerations. Needless to say that, by lowering Tr , one can easily 
reduce the predicted monopole flux by a couple of orders of magnitude below the Parker 
bound. However, this flux cannot be suppressed much further with natural (not too 
small) values of the coupling constants. In conclusion, we have shown that thermal 
inflation associated with the PQ symmetry and the \i term can naturally suppress the 
present flux of monopoles from SUSY hybrid inflation to values below but near the Parker 
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bound. This flux should be possibly accessible to ongoing and future experiments. 

The present abundance of the D and H fermions with masses ~ 1.1 x 10 12 GeV is found 
from Eq.flTTD to be Qdh ~ 0.185. These particles can therefore provide a considerable 
fraction of the cold dark matter. The rest can consist of axions which are also present in 
this scheme. The relic axion abundance fl a has been calculated in Ref. [P5fl . Assuming 
that the initial value of the axion field is about 0.55/ a , Eq.(13) of the first paper in this 
reference yields Q a ~ 0.115 for H = 65 km/s Mpc and f a = 7.57 x 10 11 GeV. We see 
that one can naturally obtain a cold dark matter component in the universe consisting 
of both axions and intermediate scale mass fermions with total energy density equal to 
about 30% of its critical density consistent with recent observations p6| , p7| . It should be 
clear that, by appropriately choosing the values of the parameters, one can easily adjust 
not only the total cold dark matter abundance but also its composition of axions and 
intermediate scale mass fermions. The present numbers only serve as an example. 

As explained, the D and H fermions are stable due to the two Z2 remnants of 
U(1)r x U(1)pq . Their stability is absolute if these t/(l)'s are exact symmetries of 
the superpotential to all orders. However, as has been shown [28], global £/(l)'s may 
be present as effective rather than exact symmetries. Indeed, some of the discrete sym- 
metries which normally emerge from the underlying string theory can effectively behave 
as continuous symmetries. These continuous symmetries are expected to be explicitly 
broken by some higher order operators in the superpotential which are allowed by the un- 
derlying discrete symmetries. If the order of these operators is adequately high, they do 
not affect our scheme except that they may provide highly suppressed Yukawa couplings 
for the decay of the .D's and/or the if's by violating either or both the Z2S. Such cou- 
plings could be DFF, DF C F C or HFF C with coefficients suppressed by (f a /mp) 5 leading 
to a lifetime ~ 10 22 years. The long-living .D's and if's with masses ~ 10 12 GeV may 
then provide an explanation [14]] of the recently observed [|I5) ultra-high energy cosmic 
ray events. Note that these intermediate scale mass particles were introduced for making 
thermal inflation possible, and thus provide a mechanism for diluting the monopoles. 
Their role in dark matter and cosmic rays is an extra bonus! 

It is generally difficult to generate the baryon asymmetry of the universe in models 
with a low reheat temperature such as our scheme. Any pre-existing baryon (or lepton) 
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asymmetry is utterly diluted by thermal inflation. Moreover, after the subsequent re- 
heating, the universe is too cold to allow baryon number violation and out-of-equilibrium 
conditions. Baryogenesis mechanisms which may be applicable here have been discussed 
in Refs. p9| , |30| . The latter uses the fact |31| that the oscillating (thermal) inflaton does 



not decay instantaneously at T r . It rather follows the usual exponential decay law. The 
"new radiation" , which is so produced, reaches a maximum temperature which is much 
higher than the electroweak scale. This radiation then gradually cools down and, finally, 
dominates the energy density at T r . During this process, a lepton asymmetry can be 
generated via the Affleck-Dine (AD) mechanism p2| . The decay of the AD condensate 
is plasma blocked at temperatures higher than its frequency of oscillations which is ex- 
pected to be of the order of the electroweak scale. Actually, this condensate decays at 
a temperature T* ~ Mw, generating a lepton asymmetry of order unity (or smaller), 
provided that its energy density at T* is comparable to (or smaller than) the "new ra- 
diation" energy density. A fraction of this asymmetry is then immediately converted 
into baryon asymmetry by the electroweak sphalerons. From T* until T r , the baryon 
asymmetry acquires |3(| a dilution factor (T r /T*) 5 and remains constant thereafter. It is 
clear that this scenario can easily lead to an adequate baryogenesis in our scheme. 

In summary, we considered the cosmological problem arising when hybrid inflation is 
applied to (SUSY) GUTs which predict magnetic monopoles. This problem is due to the 
copious monopole production at the end of inflation. We showed that the monopole flux 
can be naturally reduced to values below but near the Parker bound by invoking thermal 
inflation "driven" by the PQ symmetry, which also generates the fi term of MSSM. This 
flux may be accessible to ongoing and future experiments. Although our mechanism was 
presented within the SUSY PS model, its applicability is much wider. An interesting 
by-product is the presence of intermediate scale mass fermions and axions in the cold 
dark matter of the universe. These fermions, which were introduced for making thermal 
inflation possible, may explain the recently observed ultra-high energy cosmic rays. 

We thank A. Riotto for useful discussions. This work was supported by the Euro- 
pean Union under TMR contract number ERBFMRX-CT96-0090, the DOE under grant 
number DE-FG02-91ER40626, and the NSF under subcontract PHY-9800748. 
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